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Abstract 

We consider an open spin chain model with GL{N) bulk symmetry that is broken 
to GL(M) x GL(N — M) by the boundary, which is a generalization of a model aris- 
ing in string/gauge theory. We prove the integrability of this model by constructing 
the corresponding commuting transfer matrix. This construction uses operator-valued 
"projected" X-matrices. We solve this model for general values of N and M using the 
nested algebraic Bethe ansatz approach, despite the fact that the if-matrices are not 
diagonal. The key to obtaining this solution is an identity based on a certain factor- 
ization property of the reduced K- matrices into products of i?-matrices. Numerical 
evidence suggests that the solution is complete. 
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1 Introduction 



As shown long ago by Sklyanin pQ, the construction of an integrable open spin chain model 
requires two main ingredients: an i?-matrix (solution of the Yang-Baxter equation) which 
determines the bulk terms in the Hamiltonian, and right/left iT-matrices (solutions of the 
boundary Yang-Baxter equation j2j |3]) which determine the right /left boundary terms in the 
Hamiltonian. By now it is well understood how to solve models with diagonal i^-matrices. 
However, solving models with i^-matrices which are not diagonal in general remains a chal- 
lenging problem. Such ^-matrices can have matrix elements which are either c-numbers 
or operators. While models with non-diagonal c-number-valued i^-matrices have received 
considerable attention (see, for example, [3]-[TT]), much less is known about models with 
operator-valued .ff -matrices. 

An interesting class of non-diagonal operator- valued i^-matrices consists of so-called pro- 
jected ^-matrices found by Frahm and Slavnov [18]. Integrable open spin chains constructed 
with ^-matrices of this type have found applications in condensed matter physics [19j [20j EI] 
and string/gauge theory [221 1231 123] . R 

We consider here an integrable open spin chain model constructed with such projected 
K- matrices. The chain has L + 2 sites, labeled X, 1, . . . , L, Y. The space of states is 

XI L Y 

C M ®C N ®---C N ®C M , (1.1) 

where 1 < M < N. That is, the vector spaces of the "bulk" sites (labeled 1, . . . , L) all have 
dimension N, while the vector spaces of left and right "boundary sites" (labeled X and Y, 
respectively) have a lower dimension M. The Hamiltonian is given by 

H = Q^hjvQM + V + i + Q { y M) h L ,yQ { Y M) , (1.2) 

i=i 

where the two-site Hamiltonian hij+i is given by 

hi,i+x = li,i+x - Vi t i+i , (1.3) 

where I and V are the identity and permutation matrices on C N <g> C N , respectively; and 
Q( M ^ is a diagonal N x N matrix which projects C N to C M , 

Q M = diagQ^l^l, 1 1 ^_ I 0) . (1.4) 

M N-M 

1 For related work in string/gauge theory, see e.g. [23]-[2E] and references therein. 
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We drop the null rows and columns of the left and right boundary terms in the Hamiltonian, 
which therefore should be understood as MN x MN matrices acting on C M ® C N and 
C N <g> C M , respectively. 

Although the bulk terms have GL(N) symmetry, the boundary terms reduce the sym- 
metry to GL(M) x GL(N - M). We shall refer to this model as the GL(N)/(GL(M) x 
GL(N - M)) model. The case (N,M) = (3,2) was recently studied (following [221 123]) in 

Within the quantum inverse scattering method, the standard approach for solving in- 
tegrable spin chains with higher-rank symmetry is nested algebraic Bethe ansatz (ABA) 
[32| [33] . This approach has been adapted to open spin chains with diagonal if-matrices 
in [MJ |35j [36] . We further adapt this method, along the lines in [191 [2D1 IZE] for a related 
model with M = 2, to solve the GL(N) / (GL(M) x GL{N - M)) model for general values 
of N and M. The identity (14.3)) . which relies on a certain factorization property of the 
"reduced" i^-matrices into products of i?-matrices, plays an essential role in obtaining a 
solution. Numerical evidence suggests that the solution is complete. 

The outline of this paper is as follows. In Sec. [2] we construct the transfer matrix 
corresponding to the Hamiltonian fll.2p . thereby proving the integrability of the latter. In 
Sec. [3]we consider, as a warm-up, the special case M — 2. We establish our notation, present 
the nested ABA solution, and provide some evidence of its completeness. We then treat the 
general case M > 2 in Sec HI Finally, in Sec. [5] we present our conclusions and list some 
interesting unresolved questions. Appendix [A] contains our proof of the important identity 

2 Transfer matrix 

As already mentioned in the Introduction, the transfer matrix is constructed from an R- 
matrix and right/left K- matrices. The former is a solution R(u) of the Yang-Baxter equation 
(YBE) 

R 12 {ui - u 2 ) R\z{ux) R 2 3(u 2 ) = R 2 z{u 2 ) i?i 3 («i) R\z(ux - u 2 ) . (2.1) 

2 The idea of breaking a symmetry down to a subgroup by boundary interactions has recently been 
explored a great deal in the 0(N) case at the critical point in [29]-[3T]. 
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In view of the GL(N) symmetry of the bulk terms of the Hamiltonian, we take the well-known 
rational solution 

N N N 

R{u) =uI + iV = a(u) £ e2p ® e £) + &(„) £ e M g, e g0 + i £ e W ® e £> , (2.2) 

GL=1 a,b— 1 a, 6=1 

where 

a(w) = u + i , b(u) = u , (2.3) 

and e^P is the standard elementary N x N matrix whose (a, b) matrix element is 1, and all 
others are zero; i.e., [e^]„ = 5 ai 5 bj . 

The right i^-matrix K~(u), which here acts on C N ® C M , is a solution of the right 
boundary Yang-Baxter equation (BYBE) [U El E] 

-R12O1 - u 2 ) Ki 3 (ui) Rii(u\ + u 2 ) K^{u 2 ) 

= K^(u 2 ) R X2 {ux + u 2 ) K^M Rx 2 {u x - u 2 ) . (2.4) 

We take the solution [18] @ 

M M N M 

K-(u) = aM £ e[V ® eg° + E «2° ® 4f + «s(«) £ E e ^ ® e ^ , (2-5) 

a, 6=1 a,6=l a=M+l 6=1 

where 

ai(u) = 1 — u 2 , a2(w) = -2m , a 3 (it) = 1 + u 2 . (2.6) 

For the case (N, M) = (3, 2), this matrix coincides with the one we used earlier in [5]. (There 
we called the left and right i^-matrices K L and K R instead of K + and K~\ and we labeled 
the left and right spaces and L + 1 instead of X and K, respectively.) Unfortunately, we 
were unaware of [18] at that time. 

Here we define the left i^-matrix K + {u) to also act on C N ® C M . @ It satisfies the left 
BYBE H 

flia(-«i + u 2 ) K+iuiY 1 R 12 (- Ul -U2-V) K&M* 

= K+(u 2 ) t2 R 12 {- Ul -U2-v) KUuiY 1 Ri 2 {- Ul + u 2 ) , (2.7) 

3 See Eq. (3.12) for itiK-(u)tti in [18]. Here we take the constant c = (in order to match with the 
Hamiltonian (jl.2|l ). we rescale the rapidity u i— > — m (in order to match with our conventions for the i?-matrix 
(|2.2|) ). and we clear the denominators by performing an overall rescaling. 

4 In [S], we instead defined the left if -matrix to act on C M ® C N (with (iV,M) = (3,2)); i.e., the two 
if-matrices are related by permutation. 
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where ij denotes transposition in the i th space, and rj = iN appears in the crossing-unitarity 
relation 

R 12 (u) h R 12 (-u-ri) h ocl, (2.8) 

where the proportionality factor is some scalar function of u. A solution is provided by the 
"less obvious" isomorphism pQ 

K+(u) = tr 2 V 12 R l2 {-2u - v)K 23 {u) , (2.9) 

which gives (up to an irrelevant overall factor) 

M M N M 

ir + (t*) = 6i(«)E e 2P® e 2 +62(«)E e 2°® e £ 0+ 63(«) E E^®^ 2 - 10 ) 

a,b=l a,b=l a=M+l 6=1 

where 

b^u) = u 2 + i{N - M)u , b 2 (u) = 2iu-N, b 3 {u) = -u 2 - iMu . (2.11) 

Again, for the case (N, M) = (3, 2), this solution agrees with the one used in [5]. 
The transfer matrix t(u) is given by pQ 

t(u)=ti a K+ x (u)T al ... L {u)K- Y (u)f al ... L (u) } (2.12) 

where the trace (tr) is over an N- dimensional auxiliary space denoted by a. The argument 
of the trace acts on 

a X 1 L Y 

C N ®C M ®C N ®---C N ®C M , (2.13) 

and therefore t(u) acts on (11.11) . as does the Hamiltonian. The monodromy matrices T and 
T are given by 

T al ... L (u) = R al (u) ■ ■ ■ R aL (u) , f al ... L {u) = R aL {u) ■ ■ ■ R al {u) . (2.14) 

Indeed, it can be shown that the transfer matrix (I2.12p obeys the fundamental commutativity 
property 

[t(u),t(v)] = 0. (2.15) 
It can also be shown that this transfer matrix contains the Hamiltonian (11.21) . 

H = Cl ^-t(u) +C2I, (2.16) 



where 



C '=2iV ( - i; 



c 2 = L + 1 



N 



The relations (I2.15P - ( 12 . 1 T[) demonstrate the integrability of the Hamiltonian. 
The transfer matrix has the GL(M) x GL(N — M) symmetry 

[tH,f)®0 0L ®f)] =o, 

where 

[) G GL(M) , [)' G GL(N - M) . 




(2.17) 



(2.18) 



(2.19) 



3 Nested ABA for M = 2 

We now proceed to diagonalize the transfer matrix of the GL(N) / (GL(M) x GL(N — M)) 
model via nested ABA for the special case M = 2. 



3.1 Preliminaries 

We begin by assembling the ingredients needed to carry out the ABA analysis: suitable 
operators, pseudovacuum states and commutation relations. For the M = 2 case, the left 
X-matrix H 2 . 101) has the form (as an N x iV matrix in the auxiliary space) 

/ a n (u) a 12 {u) ■■• \ 

"2i (u) a 2 2(u) ■■• 

P(u) ■■• 

p(u) ■■■ 



K. 



aX 



U 



(3.1) 



10): 



(3.2) 



y ••• p(u) J 

where otjk{u) and P(u) are operators on the two-dimensional quantum space X. For future 
reference, we now introduce a "down" pseudovacuum state for this space, 


1 

and note that it is an eigenstate of the diagonal operators, 

an(«)|0>x = &i(«)l°>*, 
« 22 (m)|0)x = (h(u) + b 2 (u))\0) x , 
P(u)\0)x = h(u)\0) x , (3.3) 
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and is annihilated by 0:12 (u), 

a l2 (u)\0) x = 



(3.4) 



The right K- matrix (12.51) has a similar structure. Introducing a "down" pseudovacuum 
state also for the quantum space Y, 



|0>, 





1 



(3.5) 



we see that it is an eigenstate of the diagonal operators § 

[jr(u)] u |0>y = Oi(tt)|0>y, 



u 



22 



10) 



Y 



)Y 



[K-(u)]..\0) 
and is annihilated by [K~(u)} 12 , 

[K-(u 



(ai(u) + a 2 (u))\0) Y 
O3(«)|0)y, j = 3,. 



N, 



12 



\0) Y = 0. 



(3.6) 



(3.7) 



The transfer matrix (12.121) can be reexpressed as 

t{u) = tr a K+ X (u) T a l.. LY (u) 

where T~ 1 ... L y(u), defined by 



r a i 



■LY 



U 



T al ... L (u) K nY {u) T al ... L {u) 



(3i 



(3.9) 



also obeys the right BYBE (12.41) . It is from this object that we must identify suitable 
operators (among them, creation-like operators). In view of the form (13.11) of the left K- 
matrix, we follow [TU [20l [21] (see also [TQj, |34l |35l [36] and references therein) and write 
7~ al ... LY { u ) as follows (as aniVxiV matrix in the auxiliary space) 



LY 



A 



(i) 

AT-1,1 



U 



(3.10) 



5 We denote by [K (u)]j k the (j,k) element of K aY (u) considered as an N x N matrix in the auxiliary 
space, analogously to (|3.1I) . 
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where A^ k \u) , Bf\u) , Cf\u),D^\u) are operators on the quantum spaces 



1 L Y 

C N ®---C N ®c 2 . 



With respect to the all "down" pseudovacuum state 



(3.11) 



10) 



l-LY 



|0)i...i|0>. 



10) 



l-L 





V 1 / 



(3.12) 



Bf\u) and Cf ] (u) are annihilation and creation operators, respectively, 

Bf\u)\0} h .. LY = 0, Cj 1) ( M )|0) 1 ... Ly ^0, 
and .D^^m) is diagonal, 

..Ly = a(«) a 3 (M)|0)i...Ly . 



Moreover, defining the operators Aj k \u) by 



^l ) («) = Ag(u) 



a(2u) 



8 jk D^{ 



u 



(3.13) 



(3.14) 



(3.15) 



we find that 



i£>(«)|o> 1 ... Ly = 6(«) 2£ '[ir«(t»)] tt |o> 1 ...iy 



(3.16) 



where [if W (w)] are operators on the two-dimensional quantum space Y defined by 



[K-M( 



u 



[K-(u)] jk - l ^M^J. .y- A- = 1 V-l. 



a(2u) 



(3.17) 



The trace over the auxiliary space in the expression (13.81) for the transfer matrix can now 
be performed, resulting in the more explicit expression 



N-l 



t(u)=J2 [K +{1 \u)] ik A^(u) + FW(u)D^(u) 



(3.18) 



j,k=i 



where . fc are operators on the two-dimensional quantum space X defined by 



■ft 



[K + (u 



a jk (u), j,k = 1,2 
* ~ \ 6 jk p(u), j,k = 3,...,N-l 



(3.19) 



and 



a(2u) 



.3=1 



(3.20) 



Note that the expression (13.181) for t(u) does not involve either annihilation or creation 
operators, which is necessary for carrying out the nested ABA analysis. 

The operators obey the following commutation relations, 

1 



b(u — v)a{u + v) 



(u + v + i)} [J2« („_„)] ^ Cf (v) A$l («) 



a(2u)b(u — v) 



+ 



ib(2v) 



[R^(2u + i)] ijl .. k Cp\u)D^(v) 



+ 



a(2u)a(2v)a(u + v) 

a(v-u)b(v + u) (1) (1) 
b(v — u)a(v + u) J 



a(2v)b(u — v) 3 a(u + v ) J 



a^w — v) u ' jk 'j 

where R^(u) is the GL(N — 1) i?-matrix, with matrix elements 

(«)]„.,, = a(«), [««(«)].. . fc = &(«), [i? (1) ( 



(3.21) 



jk;jk 

N — 1 



(3.22) 



with a{u) and as before (12. Summation over repeated indices is understood in the 
commutation relations. 



3.2 First level 

The pseudovacuum state for the full space of states is given by 

\0)xi-ly = \0)x\0)i...ly = |0)x|0)i... L |0)y . (3.23) 

It is an eigenstate of the transfer matrix (13.181) by virtue of (13. 3p . (13.41) . (13. 6p . (13.71) . (I3.14p - 
(I3.17p . This state is not the lowest-energy state. Indeed, it is an eigenstate of the Hamiltonian 
(11.21) with energy eigenvalue 2, while there are eigenstates (such as the all "up" state) with 
energy 0. 



8 



We make the ansatz that the eigenstates of the transfer matrix (which are inde- 

pendent of the spectral parameter u by virtue of the commutativity property (12. 15ft ) can be 
obtained by acting on the pseudovacuum state with the creation operators, namely, 

= CfKi) • • K mi )|0)xi... i ^ (1),1 - i ^ , (3.24) 
where again summation over repeated indices is understood. 

By acting with the expression (I3.18P for the transfer matrix on this state, and using 
the commutation relations ( I3.21j) to repeatedly move A^'iu) and D^(u) past consecutive 
creation operators until arriving at the pseudovacuum state, two types of terms are generated. 
The "wanted" terms are those generated by the first terms in the commutation relations; 
the remaining terms are "unwanted" . The "wanted" terms give 

t(u)\n {1) ) = A(u)\n {1) ), (3.25) 

with 

mi 



AH = f (u)a(u) 2L H 



l ( u i,j — u )b{uij + u) 
11 b(u lt j — u)a{u\j + u) 



+ h ^ 2L Wu V~TZ T A (1) KK;}), (3-26) 

b[u — ui j)a[u + U\j) 



where 



f {u) = o 3 («)^63(u) + -^ T [26 1 («) + 6 2 (tt) + (JV-3)6 3 («)] 

a{2u) 



(2u + iN)(u 2 + 1] 



2 



(3.27) 



2u + i 

Moreover, A^(m ; {wij}) is a solution of the eigenvalue problem 

t^^;^!,,-}),,..^,,.,^^ 1 ^ 1 -^ ^w^Kjj/W^^, (3.28) 

where the level-one inhomogeneous transfer matrix t^'(u ; {uij}) is defined by 

t«(« ; {u ld }) = tr aW K+y x (u) T^L miY {u ; { UlJ }) (3.29) 
where now the auxiliary space, denoted by aP~\ has dimension N — 1; and 

rJl. miY (u ; {u ld }) = TX... mi ( u ; {«!„■}) ^(«) ; K,}) , (3.30) 

where the level-one inhomogeneous monodromy matrices are given by 

= ^Si) roi («-«i,™i)"--RS!) 1 («-«i,i)- (3-31) 
9 



By virtue of the fact that the level-one i^-matrices satisfy shifted BYBEs 

B${ui - u 2 ) #r 3 %i) 4? («i + U2 + i) K 23 {1 \u 2 ) 

= K^ {l \u 2 ) R^( Ul +u 2 + i) K~ (1) ( Ul ) R^( Ul - u 2 ) 



(3.32) 



R? 2 \- Ul + u 2 ) tfJ^M* 1 R^i-ux -u 2 -r)-i) <3 (1) («2)* 2 (3.33) 
= K 2 + 3 (1 W 2 - u 2 - r, - i) K+ {1) ( Ul ) h R[ 2 \- Ul + u 2 ), t] = i(N - 1) , 



(cf. Eqs. (12. 4p . (12. 7p . respectively), the level-one transfer matrix (13.291) has the commuta- 
tivity property 



[&\ui{u 1J }),t<»(v ] { Ultj })]=0. 



(3.34) 



Although for the level-one transfer matrix the auxiliary space and the "bulk" quantum spaces 
(i.e., those labeled 1, . . . , mi) have dimension one lower compared with the original transfer 
matrix, the "boundary" quantum spaces (i.e., those labeled X,Y) remain unchanged. 



3.3 Iterating 

We continue to iterate the above procedure. We define 

where the auxiliary space, denoted by aS l \ has dimension N — I, and 



>(0 



(0 



where R^(u) is the GL(N — I) .R-matrix, 

[RU(u)}.... = a(u), [R®{u)\. k . k = b{u), [R«\u)\. k]k . = i. 
j,k = l,...,N-l. 



>(0 



(0 



(3.35) 



(3.36) 



(3.37) 



We set 



(0 



j }) 



C+i), 



^AT-Z-lJV-Z-l( U ) 4v-Z-l( n ) 



|P+1) 

■i- 



l = l,...,N-3, N>4. 



(3.38) 
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The above equations are valid also for I = if we identify 

u j = , m = L , 
and also T~ {0) = T~, etc., see (Q, fl^TUj) . We define 

and find that 



(3.39) 



(3.40) 



£ (m) («)|0>i... mi y = b -^^f[a(u-u l , j )a(u + u ltj + 11)10),. 

.7=1 



■m;Y ; 



i^HlO)!...^ = Jj6(«-«ij)6(« + «w + i0 [^ (m) («)],j0)i... mi y, (3.41) 



where 



a 



ib(2u)a^(u) 



-5 jk I, j,k = l,...,N-l-l.{3A2) 



^ k b{2u + il)a{2u + il) 

The commutation relations are generalizations of (j3.2ip ; in particular, the terms which 
generate the "wanted" terms are given by 

1 



A^(u)cf^(v) 



b{u — v)a(u + v + il) 



[R^(u + v + i(l + l))] if . lh 



x[R^(u-v)] klh . kj C^ 1 \v)4i 1 \u) + 



D^(u)C? +x \v) = f- U l bi ; + U + ^ Cf +1 \v)D^\u) + .... 
1 b(v — u)a(v + u + il) J 

The level-/ transfer matrix is given by 

f w (u;{«U» = tr^K^^rJlju-Au^}) 



(3.43) 



N-l-l 



where 



and 



E [K +il+1) (u)] jk A^\u) + F^\u)D^\ 
j,k=i 



u)\ jk = [K + (u)] jk , j,k = l,...,N-l-l 



u 



(3.44) 



FV +1 \u) = /3(u) + 



a{2u + il) 



^ ajj(u) + (N — I — 3)(3(u) 



(3.45) 



(3.46) 
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The i^-matrices satisfy the shifted BYBEs 

Rikui - u 2 ) Kn {l) M + u 2 + il) K~ 3 {l \u 2 ) 

= K- {l \u 2 ) R^im + u 2 + il) K- {l \ Ul ) Rfkm - u 2 ) , (3.47) 

R? 2 (- Ul + u 2 ) K^iurf 1 i2§(-ui -u 2 -r)-il) K+ (l \u 2 ) t2 (3.48) 

= K+ (l \u 2 y* Rfl(- Ul -u 2 - V -il) K^M* Rfk-^i +u 2 ), V = i(N - > 

and therefore the level-/ transfer matrix also has the commutativity property. 
Acting with the transfer matrix (I3.44p on the Bethe state 

| (W) } = C^iu^) . . .^|( Mm , mi+1 )|0) xl ... m;+iy ^ +1 )--^ +1 , (3.49) 

the "wanted" terms give 

t®(u; K;})|^ +1) ) = A«(«; {u ld })\rt l+ V) , (3.50) 

with 

A«(« ; {«y}) = /,(«) f[ «(« - «uM« + ««< + «) II a p±M^fp±lA±l±3 

+ n6(«-ti w )6(TH-« IJ +*0 J] ^ f _„, ; ^7X^ : . n A^(«;{m + ij}),(3.51) 
i=i 

where 



^ 6(w - ui+\J)a{u + ui+ij + iZ) 



and 



fi(u) = f+(u)fr(u), (3.52) 



fr(2 M )a 3 (n) = 2u(u 2 + 1) 
/? W ~ b(2u + il) ~ 2u + il ' 

/+(«) = frfr) + * [26i(ti) + 6a(«) + (iV ~ * ~ 3) W] 

(a B + *A0(« , + D i = ,... !iV _3. (3.53) 



2w + i(Z + 1) 
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3.4 Final level 



We iterate the recursion relation (13.51 j) until we reach I = N — 3. At that stage we need the 
eigenvalues of the transfer matrix t^ -2 ) (u ; {ux-2,j}) = tr a (jv-2) K + ( N ~ 2 *> (u) T~ ( N ~ 2 > {u ; {uN-2,j}), 
where the auxiliary space a^ N ~ 2 ^ has only two dimensions. The .fT-matrices are given by 



K -(N-2) { 



U 



K +(N-2) { 



U 



( ai(u) + a 2 (u) 



V 

/ b l (u) + b 2 (u) 



\ 



ai{u) a 2 [u) 
a 2 (u) a>i(u) 



bi(u) b 2 (u) 
b 2 (u) 



ai(u) + a 2 (u) J 
\ 



6i(u) + &2(m) / 



i{N-2)a^u) 
2u + i(N - 2) 



I. 



(3.54) 



where matrix elements which are zero are left empty. They obey the shifted BYBEs (I3.47p . 
(J3HH]), respectively, with I = N -2. 

A priori, one would expect to encounter serious difficulty in diagonalizing this transfer 
matrix, since both .fT-matrices (in particular, the left one) are not diagonal. Remarkably, 
this is not the case. Indeed, we note the identity 



t^- 2 \ U] {u N ^}) 



3, 

2u 



2u + i(N 



tr a S ( aXl 2) m N _ 2 Y( U ! {^N-2,j}) S, 



{ aXl 2) m N . 2 Y( U I {uN-2j}) , (3-55) 



where 



sSi. 2) miV _ 2 y( M ; W-2,j}) 



c( N - 2 ) i r i\ 

b aX\-m N ,Y\ U \ \ U N-2, 3 )) 



= R ( a N x - 2 \u + i(N - 2)) R^- 2 \u + u N . 2jl + i(N - 2)) • • • 
x + u N -2,m N _ 2 + i(N - 2)) R^~ 2 \u + i(N - 2)) 



tl. 



(N-2) 

, aY ^^ anlN _ 2 \^ ^N- 2 ,m N _ 2 



x R ( al >(u-u N _ 2A )R i 



aX 



u 



(3.56) 



and a = a( N ~ 2 ^ is the two-dimensional auxiliary space. That is, the transfer matrix t*^ 2 ) (u ; {uj^^ 2 j}) 
is the same as the transfer matrix of an open inhomogeneous spin-1/2 GL(2)-invariant chain 
of length 2 + m N ^ 2 with trivial K- matrices (i.e., equal to the identity matrix)! A proof for 



6 A similar observation (although without proof and only for the case M = 2) has been made for related 
models in [TU [20l [2TJ. 
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general values of M is given in Appendix |A1 The corresponding eigenvalues can therefore 
be determined by standard methods such as pQ, and we obtain 

mjv-2 

A (7V ~ 2) (u ; {u N ^ 2 ,j}) = f N -2{u) n (u-u N - 2 j+i)(u + u N -2j + i(N- 1)) 

i=i 

x "jQ ( M - MiV-lj -«)(« + wjv-ij + - 2)) 



A (u - u N -xj){u + wjv-ij + ^(A 7 - 1)) 

+ /jV-l(w) ] [ I " ' • MAf-2,j)( M + M 7V-2,i + «(JV - 2)) 



J'=l 



x tt (m - u N - ltj + i)(u + un-ij + iN) , 3 

(M-U2V-lj)(u + WjV-lj+i(JV-l)) ' 



where 



/iV-2(«) 



2m(u + i) 2 (w + i(JV - l)) 2 (2u + iN) 
(2u + i(N - 2))(2u + i(N - 1)) 



2^ + *(iV-2)) 2 
= " 2u + i(N-l) • (3 - 58) 

Combining the above results, we conclude that the eigenvalues of the original transfer 
matrix (I2.12p with M = 2 and N > 3 are given by 

[ frf Qi{u + ^l) Qi+i{u + |(< + lj) 

+ fN-i{u)- — — f , (3.59) 

Qat_i(u + ^(jV - 1)) J 

where 

Qi(«) = II( w_w w)( w + w w)' i = l,...,^-l, (3-60) 

and we have made the shifts uij i— >■ My — |Z. We recall that the functions /;(«) are given by 
(O, ([333]), 0338]). 

We have thus far ignored all the contributions from "unwanted" terms in the commutation 
relations. Such contributions vanish provided the parameters {uij} satisfy the Bethe ansatz 
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equations 



ei(Mi,fc) 2L = 6i (u 1)k ) JJe 2 («i,jfc -«i j ) e 2 (wi, fc + 



/T( 2 



x JJ e_ 1 (wi )fc - w 2 j) e_i(w 1)fc + M 2j ) , k = 1 , . . . , m 1 



mi 



mi-i 



i=i 



3=1 



X 



e_i(wj j& - iij+ij) e_i(«; )fc + , A; = l,...,m z , / = 2, . . . , N - 2 , 



m N _ 1 



1 = 6jv-i(«Ar-i,fc) _Q e 2 (MAr_i )ifc - un-ij) e 2 {u N -x,k + un-ij) 



3 = 1 



mjv-2 



X I I e -l( u iV-l,fc ~~ U N-2,j) e_i(uAr_i ifc + UN-2,j) , fc = 1 , . . . , mjy-1 , (3.61) 



where 



6,M = I 



eiV-3(") 2 
ejv_i(n) 2 



l = N -2 



l = N -I 



otherwise 



(3.62) 



and we have used the standard notation 



e n \u 



u + in/2 
u — in/2 



(3.63) 



Finally, from the relation ( 12 . 16[) between the transfer matrix and the Hamiltonian, we 
find that the energy eigenvalues are given by 



d mi i 
E = c 1 —A{u) +c 2 = 2 + V^ — 



(3.64) 



3.5 The case N = 3,M = 2 

For the case (N,M) = (3,2), the above results do not coincide with those in our previous 
work [21] . Indeed, there we found that the eigenvalues are given by the same expression 



15 



(I3.59P but with different functions fi(u), namely, 



previous 
! 



U 



ff rcvious (u) 



f: 



previous , 



U 



(2u + 3i)(u + if 
2u + i 

u 3 (2u + 3i)(u + i) 
2u + i 

-u 3 (u + i) = f 2 (u). 



u + i 



u — i 



u 



fo(u) 



u + 2i 



fi(u) 



(3.65) 



(See Eqs. (2.33) and (2.36) in [21].) Equivalently, the two sets of results can instead be 
related by 



* <• -2/ % \-2 



Qi(u) = g(u)Qr vlous (u), g(u) = (u + -)~\u - -) 



Q 2 (u) = Q 



previous , 



U 



(3.66) 



since 



(•previous / 
JO V 



U 



fo(u) 



g(u+ 3 i) _ /r vious ( 



u 



g(u + 



h 



u 



(3.67) 



The discrepancy in the two sets of results arises from different choices of pseudovacua. 
In [23] we chose the pseudovacuum to be a ground state (E = 0), while here we have taken 
the pseudovacuum to be an excited state (E = 2). (Notice the additive constant in the 
expression ( 13. 64ft for the energy.) 

We have performed a numerical analysis of completeness of the new solution for small 
values of L along the lines discussed in Appendix B of |24j. The results for the case L = 3, for 
which case there are M 2 N L = 108 states, are displayed in Table [H Although we find some 
levels for which m 2 > nil (which we did not find with our previous solution), this solution 
also appears to be complete, at least for small values of L. Note that the Bethe roots for the 
ground (E = 0) state have a rather complicated structure. Comparing this table with Table 
2 in Ref. [23], we see little apparent relation between the two sets of Bethe roots describing 
a given energy level. 

It would be interesting to re-derive our previous solution [23] (obtained by analytic Bethe 
ansatz, which is a heuristic approach) by the more rigorous nested ABA approach consid- 
ered here. Unfortunately, we have so far not succeeded. Indeed, if we try to use the all 
"up" state as the pseudovacuum, then the creation and annihilation operators seem to be 



.4 



12 1 



U 



B^\u) and A^Wu) , C\ L> (u), respectively; hence the transfer matrix seems to involve 



creation and annihilation operators. 
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3.6 The cases N > 3 , M = 2 



For iV > 3 , M = 2, the solution also seems to be complete. For example, we display in Table 
[2] our results for (JV, M) = (4, 2) and L = 2, for which case there are M 2 N L = 64 states. 



4 Nested ABA for general values of M 



For M > 2, the left if- matrix has the form (cf. Eq. (13.11) ) 

/ a u (u) ■■■ a 1M (u) 



a M1 (u) 




%m(«) 
0(u) 



o \ 






(4.1) 



\ ••• ••■ /3(u) / 

We can therefore iterate the recursion relation (13.511) until we reach I = JV — M — 1. At 
that stage we need the eigenvalues of the transfer matrix t^ N ~ M \u ; {wjv-Mj})) f° r which the 
auxiliary space aS N ~ M "> has dimension M. The corresponding if -matrices are given by 

i(N -M)a 3 (u) 



K-( N ~ M \u) 

K +(N-M) (u) 



ai(u) - 



l + a 2 {u)V 



2u + i(N - M)_ 
hiujl + hiujV , (4.2) 
where I and V are the identity and permutation matrices on C M (g> C M , respectively. They 



obey the shifted BYBEs fl3TTj) . (13T48D . respectively, with l = N - M. 

Since these if-matrices (in particular, the left one) are not diagonal, it is not evident 
how to diagonalize the transfer matrix. Fortunately, there is an identity generalizing (I3.55p . 
(I3.56p . namely 



2u 



2u + i(N -M) 



tr a S ixi- M 1 N - M Y( U '> { U N-M,j}) S ( 



£i.^L_j,y( u I {un-mj}) ,(4.3) 



where 



C. M y^ ; W-mj}) = R ( a N x - M \u + z(JV - M)) R%- M) (u + ujv— m,i + z(JV - M)) • • • 

x + ^_M, mjv _ M + i(JV - M)) R%- M \u + z(JV - M)) 



5,(JV-M) 



(u ; {ujv-Mj}) 



, m N-M l 



X «ol [U ~ U N -M,l) H-aX 1 



«) 



(4.4) 
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and a = a^- M ) is the M-dimensional auxiliary space. That is, the transfer matrix is the 
same as that of an open inhomogeneous GL(M)-invariant chain of length 2 + tun-m with 
trivial i^-matrices and spins in the vector (M-dimensional) representation. See Appendix 
lAl for a proof. The corresponding eigenvalues can be found by the "ordinary" nested ABA 
[3U [35], and we obtain 

rriff-M 

A( N - M \ U ; {u N - M j}) = f N - M {u) Y[ {u~ un-mj + i)(u + U N -MJ + l{N - M + 1)) 

3=1 

niJV-M+l 



n 



u - u N - M+ ij -i)(u + u N - M +i,j + i(N-M)) 
L (u- u N - M +i,j)(u + u N - M +i,j + i{N - M + 1)) 



mjv-M 



+ I| (w - u N - M ,j){u + u N - M ,j + i(N - M)) 

3=1 



N-2 



nil 



(u - u\j +i)(u + uij + i(l + 1)) "yj (u - ui+ij — i)(u + w/+ij + il) 



n 



Z=AT-M+1 



(w - + «/j + il) (w - U| + i,j)(tt + + i(Z + 1)) 



(u - Un-ij +i)(u + Mjv-ij + i-^V) 



n (u-ujf.^iu + us-u + KN-l)) 



(4.5) 



where 



2m(m+») 2 (u+i(N-M+l)) 2 (2u+iN) j _ at _ i\ r 
(2u+i(N-M))(2u+i(N-M+l)) 1 JV m 



(4.6) 



2M 3 (M+i(V-Af)) 2 (2M+tV) 
(2w+a)(2u+i(J+l)) 



I = N-M + 1,...,A- 1 



Combining this result with those from the recursion relation (I3.5ip . we conclude that the 
eigenvalues of the original transfer matrix (12.121) are given by 

N-2 



A(«) 



= fo(u)(u + i) 
+ /jv-i(w 



2L Qi(u 2 ) ^ 2L 

g^i(^ + l(A + i)) 

Q^i^ + KA-l)) 



q,(» + 1 (j + 2)) Qi + i(u + ±(i-i)) 

Qi(u + U) Q l+l {u + Ul + l)) 



(4.7) 



where 



( 2u(u 2 + 1) 2 (2u+iN) 
(2u+il)(2u+i(l+l)) 



l = 0, 



N-M-l 



2u(u+i) 2 (u+i(N-M+l)) 2 (2u+iN) i _ at _ M 
(2u+i(N-M))(2u+i(N-M+l)) 1 ~ JV 1V1 



2u A (u+i(N-M)) 2 (2u+iN) 
(2u+il)(2u+i(l+l)) 



I = N - M + 1,...,N -1 



Ill, 



Qi{ u ) = - uij)(u + uij) , I = 1, . . . ,N — 1 
i=i 

and (as before) we have made the shifts uij uij — ~£. 
The corresponding Bethe ansatz equations are given by 

mi 

eiOi,fc) 2L = ®i(ui,k)\\e 2 (u 1)k -uij)e2{u 1 ,k + u 1)j ) 



3=1 



(4.9) 



x JJ e_i(wi )fe - w 2 ,i) e_i(«i )Jfc + m 2 . 



fc = 1, . . . ,171!, 



3=1 



rn, 



mi-i 



mi+l 



X 



Y[ e-i{ui t k - ui+ij) e-.^uw + ui+ij) , k = l,...,mi, l = 2,...,N-2, 



mjv_i 



1 = 9jv-i(wjv-i,fc) e 2 («jv_i,fc - «jy-i,i) e 2 (wAr-i,fc + Wjv-ij) 



3=1 



TtiV-2 



X IT e -l( U ^-l,fe - Un-2,j) e_i(MAr_i ifc + U N ~2,j) , fc = 1 , . . . , 771JV-1 , (4.10) 



where now 



' ej v_A/+ 2 (w) 2 l = N-M 



ejy-Af-i(M) 2 
e]v-M+i(M) 2 



/ = N-M+ 1 



otherwise 



(4.11) 



and e n (w) is defined in f 1 3 . 6 3 j) . The energy eigenvalues are given by the same formula (I3.64p . 
The identity ( 14. 3ft . the expression (14. 7j) - ( j479|) for the eigenvalues of the transfer matrix and 
the corresponding Bethe ansatz equations (14. lOj) . (14.111) are the main results of this paper. 

For M > 2, this solution also seems to be complete, as is the case for M = 2 discussed 
in Sec. 13.61 For example, we display in Table [3] our results for (N, M) = (4, 3) and L = 2, 
for which case there are M 2 N L = 144 states. 
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5 Conclusions 



We have considered the GL(N)/(GL(M) x GL(N - M)) model with Hamiltonian fOD . 
which is a generalization of a model arising in string/gauge theory. We have proved the 
integrability of this model by constructing the corresponding commuting transfer matrix. 
The latter makes use of the non-diagonal operator- valued K- matrices found in |18j . 

We have found a Bethe ansatz solution of this model for general values of N and M 
using the nested ABA approach, despite the fact that the i^-matrices are not diagonal. 
The main results are the eigenvalues (I4.7p - fl4.9p and Bethe ansatz equations f l4.10p . (14. lip . 
The key to obtaining this solution is the identity (14. 3p . which relies on the factorization 
property flA.ip of the "reduced" (level N — M) i^-matrices into products of .R-matrices. In 
hindsight, this property is not too surprising, since the projected i^-matrices originate from 
"dressed" diagonal if-matrices [H]. For the case (N,M) = (3,2), this solution is not the 
same as the one found in [21] using analytic Bethe ansatz, as the two solutions are based on 
different pseudovacua. Nevertheless, numerical evidence suggests that both N = 3 solutions 
are complete. Moreover, the nested ABA solution appears to be complete for general values 
of N and M. 

Many interesting questions remain unanswered. It is unusual for an integrable model 
with a non-graded symmetry algebra to have more than one Bethe ansatz solution. (Models 
with graded symmetry algebras are known to have more than one Bethe ansatz solution, 
corresponding to the non-uniqueness of the associated Dynkin diagrams. See e.g. [3 TJ and 
references therein.) This underscores the question of whether the two proposed solutions for 
the case (N, M) = (3, 2) (namely, the one found in [21] by analytic Bethe ansatz, and the one 
found here by nested ABA) are equivalent. As noted in Sec. I3.5[ one would like to have a 
more rigorous derivation of the solution found in [24J . Similarly, for general values of N and 
M, there may be additional equivalent solutions based on different pseudovacua. Perhaps 
Bethe ansatz equations for generic open spin chains (or at least for open chains constructed 
with projected i^-matrices) can be formulated in terms of group theory data (namely, the 
"bulk" symmetry algebra and the unbroken "boundary" symmetry subalgebra); and the 
multiplicity of Bethe ansatz solutions reflects the various ways of choosing the boundary 
symmetry subalgebra. We hope to be able to address these and related questions in the 
future. 
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A Proof of the transfer-matrix identity 

Our proof of the transfer-matrix identity (13.551) . H4.3[) is based on the following remarkable 
factorization property of the "reduced" i^-matrices (i.e., the (N,M) projected K- matrices 
at level N — M) into products of .R-matrices, 

K: r- M \u) = - - - , ( 2 ; _ u) R'Sr M \u + ,(n - m)) fl g- M) M . 

K:r- M \u) = tnP^R%r M \u + i(N-M))Rlr M \u), (A.l) 

which can be verified from the expressions (13.541) . (14. 2p . Omitting the quantum-space indices 
and denoting the M-dimensional auxiliary space by a in order to streamline the notation, 
we have 

t ( N - M \u) = ti a Kt {N - M \u)T- {N - M \u) 

= tr ab V ab J#™> (« + i(N - M)) R { b N - M) (it) T~ (u) 
= tr ab V ab itf - M \u + i(N - M)) T~ iN - M \u) itf ~ M) («) 
= tr a Rf ~ M \u + i(N - M)) r~ {N - M \u) R { a N - M) (u) 

= tr a R^- M \u + i(N - M)) T^ N - M \u) K;^ N - M \u) ff~ M \u) R^- M \u) 

x Ri N ~ M \u + i{N - M)) R^- M \u) ff~ M \u) Rf- M \u) 

= ^ rtr a Si N - M \u)Si N - M \u), (A.2) 

2u + i(N-M) a K ' a v ; ' K ' 

where in the last line we have used the fact (see Eqs. (14. 4p . (13.361) ) 

Si N ' M \u) = Rl N - M \u + i(N -M))T^ N - M \u)R^- M \u + i{N - M)) , 
S[ N ~ M \u) = R[ N ' M \u)fj 1 N - M \u)R^- M \u). (A.3) 
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E 


s 


K*} 







5/2 


1. 11803?', 0.442686 ± 1.0936* 




0.381966 


3/2 


1.2944i, 0.375279 ± 1.36374? 





0.585786 


2 


0.204205 ± 1.22426? 


— 


0.82259 


1/2, 1 


0.15313 ± 1.36461? 





1.07919 





1.36676?, 1.88488? 


0,1.56857?' 


1.26795 


3/2 


1.27123? 


— 


1.38197 


3/2 


0.936268 , 0.180565 ± 1.20371? 





1.38197 


1/2 


1.36676? 





1.58579 


1/2, 1 


1.91214, 1.31987?' 





1.69722 


1/2 


1.88488?' 





2 


2 


0.866025,1. 11803?' 


— 


2 


1,3/2 


— 


— 


2 


0, 1 


— 





2.58579 


1 


0.639467,1.15027?' 


1.0505 


2.61803 


3/2 


0.322878 ± 0.500421?' , 1.04607? 





3 





0.606658,1.36676?' 


0,0.707107?' 


3 


1/2 


0.866025 


1 


3.31526 





0.606658,1.88488?' 


1.15861? 


3.32164 


1/2, 1 


0.451092,1.17552? 





3.41421 


2 


0.479032 ± 0.521886?' 


— 


3.61803 


1/2 


0.606658 





3.61803 


3/2 


0.331608 , 0.404442 ± 0.90768?' 





4 


1 


±0.5?' 


1 


4.41421 


1/2, 1 


0.0774471 , 0.959277?' 





4.68474 





0.229729,1.36676?' 


0,0.810943?' 


4.73205 


3/2 


0.340625 


— 


5 


1/2 


0.288675 


0.745356 


5 





0.229729,1.88488?' 


0,1.22474?' 


5.30278 


1/2 


0.229729 





5.41421 


1 


0.301797,1.35023 


0.62964 


5.85577 


1/2, 1 


0.248411,1.13757 





6.92081 





0.229729,0.606658 


0,0.678531 



Table 1: Energy, spin, and Bethe roots for N = 3 , M = 2 , L = 3. 
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E 


deg 


{ui,k} 


{U2,k} 


{ U 3,k} 





5 


0.238862 ± 0.9867732 


0.240994 ± 1.546422 


- 


0.585786 


3 


0.204205 ± 1.224262 


±1.847762 


1.398972 


1 


8 


1.11803z 


1.63299i 


- 


1.26795 


5 


1.27123i 


1.797792 


1.019152 


2 


9 


- 


- 


- 


2 


7 


- 


- 


0.866025i 


2 


3 


0.866025, 1.11803z 


0.8556,1.65289* 


0.866025i 


3 


8 


0.866025 







3.41421 


3 


0.479032 ±0.521886i 


±0.765367i 


0.736813i 


4 


7 


0.5 






4 


1 


0.5 




0.866025i 


4.73205 


5 


0.340625 


0.481717 


0.6792092 



Table 2: Energy, degeneracy, and Bethe roots for N = 4 , M = 2 , L = 2. 



E 


deg 


{Ul,k} 


i U 2,k} 


{«3,fc} 





15 


0.238862 ± 0.9867732 






0.585786 


15 


0.204205 ± 1.224262 







1 


10 


1.118032 






1.26795 


14 


1.271232 







2 


6 








2 


11 









2 


15 


0.866025,1.118032 







2.58579 


3 


0.639467,1.150272 


0,1.0505 


0.89542 


3 


10 


0.866025 






3 


1 


0.866025 


0, 1 


0.866025 


3.41421 


15 


0.479032 ± 0.521886i 







4 


8 


0.5 


0.816497 




4 


3 


±0.52 


0,1 


0.866025 


4.73205 


14 


0.340625 







5 


1 


0.288675 


0, 0.745356 


0.726483 


5.41421 


3 


0.301797,1.35023 


0,0.62964 


0.669495 



Table 3: Energy, degeneracy, and Bethe roots foriV = 4,M = 3,L = 2. 
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